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Abstract 


A detailed analysis and examples are presented that show how to enrich the kinematics of 
classical Kirchhojf plate theory by appending them with a set of continuous piecewise-cubic 
functions. This analysis is used to obtain functions that contain the effects of laminate 
heterogeneity and asymmetry on the variations of the inplane displacements and transverse 
shearing stresses, for use with a {3, 0} plate theory in which these distributions are specified 
apriori. The functions used for the enrichment are based on the improved zig-zag plate theory 
presented recently by Tessler, Di Scuva, and Gherlone. With the approach presented herein, the 
inplane displacements are represented by a set of continuous piecewise-cubic functions, and the 
transverse shearing stresses and strains are represented by a set of piecewise-quadratic functions 
that are discontinuous at the ply interfaces. 


Introduction 

Many variants of refined theories for laminated-composite and sandwich plates have 
appeared in the technical literature for many years . 110 Each of these theories generally have a 
different degree of complexity and a corresponding range of validity. Plates that exhibit a 
relatively large amount of transverse shearing flexibility typically require a refined theory in order 
to obtain accurate predictions of their structural response, particularly when the characteristic 
dimension of the response is on the order of the plate thickness. In contrast, for plate-response 
phenomena that consist of deformations with relatively large characteristic dimensions, simpler 
theories that are less complex and computationally expensive can be utilized to obtain useful 
results. 

The most basic plate theory that accounts for transverse shearing deformations is known as 
the first-order shear-deformation theory (FSDT). With regards to its origins, it is also commonly 
referred to as the Reissner-Mindlin plate theory. 11 " 1 ’ In this theory, the deformation of a plate is 
expressed in terms of three translational displacements and two rotations associated with through- 
the-thickness shearing, in contast to the classical plate theory (CLPT) conceived by Kirchhoff that 
uses only the three translational displacements. A drawback associated with the relative simplicity 
of FSDT is that the through-the-thickness distributions of the transverse shearing strains are 
approximated as uniform distributions. These approximate strain distributions correspond to 
transverse shearing stresses that violate the traction boundary conditions on the top and bottom 
bounding surfaces of a plate. As a result, a shear-correction factor must be used to obtain adequate 
results. 

One step up from FSDT, is a group of relatively popular plate theories that satisfy the traction 
boundary conditions on the top and bottom bounding surfaces of a plate and avoid the need for a 
shear-correction factor . 14 41 This boundary-condition problem is mitigated in these theories, at the 
expense of more complexity, by appending the kinematics of CFPT with terms that account for 
transverse shearing by using two additional kinemetic variables and presumed distributions of 
through-the thickness shearing strains. Typically, these distributions are presumed to be 
symmetric over the plate thickness, even for laminated plates with asymmetric ply distributions. 
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Although, variations in actual shape of the transverse-shearing strain distributions for moderately 
thick laminated and sandwich plates may have a relatively small effect on the prediction of global 
response phenomena like buckling loads and fundamental frequencies, the mismatch between the 
distribution symmetry and ply lay-up symmetry is a fundamental issue that should be addressed. 
This issue is addressed in the present paper by enriching the CLPT kinematics with zig-zag 
variations to obtained transverse shearing strain distributions that are based on the actual through- 
the-thickness construction of a plate and without introducing any aditional kinematic variables. 
To accomplish this task, the baseline plate kinematics are presented first. Then, the details of the 
procedure used to enrich the kinematics are given. Lastly, examples are presented for an 
unsymmetric laminate and a sandwich plate with an isotropic core and two identical laminated 
face plates. 


The Baseline Plate Kinematics 

Consider a uniform-thickness plate and the coordinate system shown in figure 1. Points of the 
plate are located by the orthogonal Cartesian coordinates (x, y, z), where (x, y, 0), 0 < x < a, and 
0 < y < b define points of the rectangular midplane. The plate lengths in the x- and y-directions 
are denoted by a and b, respectively. Points above and below the midplane are given by the 
nonzero z-coordinate values within -h/2 < z < h/2, where h is the plate thickness dimension. The 
principal material coordinate system at the point (x, y, z) is shown in figure 2. For laminated 
plates, the fiber orientation is denoted by the angle 0(z) shown in figure 2. 

The baseline plate kinematics used in the present study are defined by 


U(x, y, z) 

= u(x, y) Z 3x + A x (z) H xz (x, y) 

(la) 

V(x, y, z) 

= v ( x ’ y ) z dy + A y (z) H yz (x,y) 

(lb) 


3 

* 

N 

II 

* 

(1c) 


where U(x, y, z), V(x, y, z), and W(x, y, z) are the displacement-field components in the x-, y-, 
and z-coordinate directions, respectively, of the material point (x, y, z). The functions u(x, y) and 
v(x, y) are the corresponding inplane displacements of the material point (x, y, 0) of the midplane, 
and w(x, y) is the out-of-plane displacements of the material point (x, y, 0). The functions A x (z) 
and A y (z) are selected to satisfy the traction-free boundary conditions on the transverse shear 

stresses at the bounding surfaces of the plate given by the coordinates (x, y, ± b) . In addition, 
A x (z) and A y (z) are required to satisfy the conditions U(x, y, 0) = u(x, y) and V(x, y, 0) = v(x, y), 
but are otherwise arbitrary. The unknown functions H xz (x, y) and H yz (x, y) account for the effects 
of transverse shear deformation. Inspection of equations (1) reveals that CPT is contained as a 
proper, well-defined subset that is obtained by setting A x (z) and A y (z) equal to zero. Several 
choices for A x (z) and A y (z) have appeared in the technical literature; for example, see reference 
32. When A x (z) and A y (z) are specified as cubic polynomials in z, the plate theory is sometimes 
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referred to as a {3, 0} theory, which indicates that the inplane displacements U(x, y, z) and 
V(x, y, z) are cubic and the out-of-plane displacement W(x, y, z) is constant across the plate 
thickness. 

The linear strain-displacement relations used in the present study are obtained by substituting 
equations (1) into linear strain-displacement relations of the linear theory of elasticity, given by 


exx( W )=f 

(2a) 

u^y-)=f 

(2b) 


(2c) 

Vj*- **)=% + % 

(2d) 


(2e) 


(2f) 


First, consider the transverse shearing strains defined by equations (2e) and (2f). Substituting the 
presumed displacement field, given by equations (1), into these strain-displacement expressions 
gives 

/ Y»(w) \ = f A,'(z) H„(x, ,) \ 

\y „( x -y- z ) I \ A,'(z) HJ.'.y) j 

At the plate midplane z = 0, the transverse shearing strains are given by 

/ Y°z( x ’y) \ = / Yxz( x >y>°) \ _ f A,'(°) H sz (x, y) j 
\Yyz(x^y)/ \ Yy Z ( x ’ y’ °) / | A/(0)H yz ( X ,y) J 

Using equation (4), the inplane displacement field is expressed as 

U(x, y, z) = u(x, y) - z aW g x ^ + F x (z) y° z (x, y) 

V(x, y, z) = v(x, y) - z 0W gy + F y (z) y° z (x, y) 

where 


(3) 

(4) 

(5a) 

(5b) 
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(5c) 


F x (z) = 


A x (z) 

A x '(0) 


F y (z) = 


A^Z) 

A/(0) 


(5d) 


In these equations, A/(0) * 0 and A/(0) * 0 are required, and F x (0) = F y (0) = 0 are required 

to obtain U(x, y, 0) = u(x, y) and V(x, y, 0) = v(x, y). Substituting equations (lc), (5a), and (5b) into 
the remaining strain-displacement expressions yields the following expressions for the nonzero 
strains 


( e xx (x, y, z) 

l ^yy( X ’ Y’ Z ) 
l Yxy( X >y>z) 


and 


du(x, y) 


f 2 \ 

a W(x, y) 

dx 


dx 2 

) dv(x,y) 

_ J 

d w(x, y) \ 

\ dy 


a 2 

dy 

dv(x, y) du(x, y) 


0 d"w(x, y) 

y dx + dy j 


V 3x ^y j 


F x (z) 

F y (z) 

F y (z) 


dyl(x, y) 
dx 

K 

dy 

dyl( x ’ y) 

dx 


y) 


+ F x (z) 



( 6 ) 


/ Yxz(x,y,z) \ = / F/(z) Y : z (x,y) \ 
\ Y yz (x, y, Z ) / |F y '(z) y ;(x,y) | 


Next, to get the complete picture of the plate deformation, as shown in figures 3-7, consider 
the components of the linear rotation vector given by 



Substituting equations (lc), (5a), and (5b) into these expressions yields 


I > dw(x, y ) i,o, > 

“x( x >y,z) = — ^ yF y (z) Y y2 (x, y) 


2 y 


dw(x, y) 


CD y (x, y, z) = ^F x '(z) Yxz (x, y) dx 


C0 2 (x, y, z) = T 


av(x, y ) 3u(x, y ) i a Y yz( x -y) c 3 yL( x ’ y) 

*y\ Z / *x\ Z / 


dx 


dy 


dy 


(8a) 

(8b) 

(8c) 


(9a) 

(9b) 

(9c) 


4 



At the plate midplane, 


I n \ \ 3w (x,y) I 

0 ) x (x, y, 0) = tO x (x, y) = ^ ^yjx, y) 

( n \ o( \ 1 o / \ <?w(x, y ) 

CD y (x, y, 0) ^ CO y (x, y) = ^y xz (x, y) ^ 


(O z (x, y, 0) =tt» z (x, y) = 2 


| f dv(x, y) du(x, y) 


dx 


dy 


(10a) 

(10b) 

(10c) 


where F x (0) = F y (0) = 0 have been used and, as shown by equations (5c) and (5d), it is noted that 
F/(0) = F/(0) = 1 . 


The constitutive equations used in the present study are those for a plate made of one or more 
layers of linear elastic, specially orthotropic materials that are in a state of plane stress. These 
equations, referred to the plate (x, y, z) coordinate system are given by 


Oxx( x > z ) 
a yy (X, y. Z) 
Cfxy( X > y> z) 


_(k) _(k) _(k) 

Qn Q, Q,6 

— (k) — (k) — (k) 

Q,: Q 22 Q,6 

— (k) _(k) — (k) 

Q 16 Q 26 Q 66 


£xx( x > y> z ) 
£ yy( x ’ y> z ) 

Yxy( X > y^z) 


(11a) 


/ C'yz^y^z) \ 
\ a xz ( x ,y,z) j 


( k) ( k) 

^44 ^45 
-ik) — (k) 


/ Y yz( X ’ y. z) \ 

\ Yxz( x ^ y^ z) / 


(Hb) 


where the superscript (k) denotes the k th ply of a laminated plate with a total of N plies. The 

— (k) 

Qij terms are the transformed, reduced (plane stress) stiffnesses of classical laminated-plate 
theory, for the k ply, and the C ;j terms are the corresponding stiffnesses of a generally 

orthotropic ply. For a heterogeneous, laminated plate with N plies, both the Qy and the C ,j 

terms are functions of the through-the-thickness coordinate z, as indicated by the superscript (k). 
In particular, the fiber orientation of the kth ply within a laminated plate is given by the angle 0(z), 
shown in figure 2, where the values of z correspond to points within the thickness of the k th ply. 
From equation (1 lb), it follows that the transverse shearing stresses on the top and bottom plate 
surfaces vanish provided that the corresponding transverse shearing strains vanish. Thus, from 
equation (7), it also follows that 


F/ 




= 0 


( 12 ) 


is required for the transverse shearing stresses on the top and bottom plate surfaces to vanish. 


5 



Equations (lc), (5a), and (5b); the conditions defined by equation (12); and F x (0) = F y (0) = 0 
define completely the kinemetics of the baseline theory used herein. The kinematics of CFPT are 
recovered from these equations by specifying F x (z) = F y (z) = 0. Fikewise, the kinematics of FSDT 
are recovered from these equations by specifying F x (z) = F y (z) = z such that equations (5 a) and 
(5b) become 


U(x, y, z) = u(x, y) + 


yl( x ’y) - 


3w(x, y) 


<3x 


(13a) 


V(x, y, z) = v(x, y) + z 


• y) 


()w(x, y) 

dy 


(13b) 


The corresponding FSDT kinematics, in the x-z plane, are illustrated in figures 3 and 4. In 
contrast, the kinematics in the x-z plane of the baseline theory presented herein are illustrated in 
figures 5-7. 


Zig-Zag Enrichment of the Kinematics 

In the baseline kinematics, the functions F x (z) and F y (z) are typically specified as 
continuous functions, with continuous derivative, such as those given in reference 32. Moreover, 
for the most part, the functions are specified as antisymmetric functions; that is, F x (- z) = - F x (z) 
and F y (-z) = - F y (z). This antisymmetry with respect to the plate midplane appears to be 
inconsistent with laminates that exhibit a high degree of asymmetry in their ply constitutive 
properties. To obtain expressions for F x (z) and F y (z) that are based on actual laminate 
construction, particularly for asymmetric laminates, these functions are partitioned herein as 
follows 


F x (z) = H x (z) + o‘ k) (z) (14a) 

F y (z) = H v (z) + <I> y k| (z) (14b) 

where H x (z) and H y (z) are continuous functions, with continuous derivatives, such as those given 
in reference 32 which include linear, cubic, trigonometric, hyperbolic, and exponential functions. 

The terms of’fz) and <J> y ’(z) are piecewise linear functions with piecewise constant 
derivatives that are referred to herein as zig-zag functions. In these equations, the superscript (k) 
also denotes the ply number of a laminate composed of N plies. Moreover, for the k th ply, it 
follows that z (kI) <z< z (k) , where z (kll and z (k) are the ply coordinates shown in figure 8 for the 
x-z plane. Substituting equations (14) into equations (5a) and (5b) yields the inplane 
displacements 


U(x, y, z) = u(x, y) — Z 


dw(x, y) 
dx 



(15a) 
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V(x, y, z) = v(x, y) - z dW j*’ ^ + [H y (z) + <’(z)] Y ;(x, y) 


(15b) 


The user-specified functions H x (z) and H y (z) define a smooth distribution of displacements 

associated with transverse shearing deformations and of ’(z) and 'If (z) are corrections 
associated with laminate heterogeneity and asymmetry. 

Following the refined zig-zag theory of Tessler et. al., 43 ply interface displacements A/ 

A W 

(see figure 8) and A y are defined in terms of the zig-zag functions by 


, (o 





(16a) 

(16b) 


for k = 0, 1, N; with aJ; 1 = A y 1 = a! 1 = A y 1 = 0 . Within the k th ply, z (k _ n < z < z (k) and linear 
interpolation gives 


where 


«(k), , .(k-l) „(k)/ \ 

Ox (z) = A x + |3 X (z - z,^)] 

• W, x A ( k -U Q (k)/ \ 

O y (z)=A y + by (z - Z (k _j , ) 



A ,k) A ( k -‘) 

A - A 

x x 



A |k| A |k " 1 ) 

A y -A y 


(17a) 

(17b) 

(18a) 

(18b) 


and where h (k) = z (k) - z (k l) is thickness of the k th ply, as shown in figure 8. For z (k l) < z < z (k) , it 
follows that the piecewise-constant derivatives of the zig-zag functions are given by 

^'(z) = pf 

or(z)= P r 

Thus, differentiating equations (14) gives 

F x '(z) = Hf(z) + pf’ ( 20a) 


(19a) 

(19b) 
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(20b) 


F/(z)=H/(z) + p 


Substituting these results into equation (7) yields 


/ Yxz(*,y,z) \ = 
\Y y z(x,y,z)/ 


H x '(z) + P x ' Y° z ( x >y) | 
H /( z ) + Py k| Yyz( x >y) f 


( 21 ) 


for the transverse shearing strains. 


At this point in the analysis, it is convenient to use equation (21) to express constitutive 
equations (1 lb) as 


/ °yz( x > y> z )\ 
\a xz (x, y,z) | 


pW p |k) 
^44 ^45 

'-'AS SS 


H y '(z) + Pf 

Yyz( X -y)\ 

H/(z) + p ( : r 

Yl( x -y) j 


(22) 


for laminated plates, where C 44 , C 45 , and C 55 are constant-valued for material points of the k th 
ply located by z (k l) < z < z (k) . To proceed following the refined theory of Tessler et. al., 43 the 
transverse shearing stresses are partitioned into 


where 


/ O yz (x, y,z)\ 
\o xz (x, y, z) f 


— (k) — (k) 

c c 

^44 v -'45 

— (k) — (k) 

c c 

^45 ^55 


H,'(z) - l]Y«(x.y)( / a"d 

h/(z) - l y"z( x> y} J \ a -/ 


K \ 

\ a -l 


TV 1 

''-'44 ^45 

pM r {k] 

SS 



Next, the ply parameters 



(23a) 


(23b) 


(24a) 

(24b) 


(k) (k) (k) (k) 

are introduced to eliminate P x and (> y from equation (23b). Because P x and P y appear as 

(k) 

independent unknowns in the analysis, they can be normalized by requiring G x = G x and 
Gy’ = G y for k=l,2, ..., N such that 
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(25a) 


R (k, = -^-l 

r* (k) 

c 


and 


c ,k| 


« G y i 

y pW 

^44 

(25b) 

-( k ) 


( x ’ y) + G ,z|Y>y) 

(26a) 

c 44 


y Yxz(x,y) +G y Y° 2 (x,y) 

(26b) 


Following Tessler et. al., the constants G x and G v are found by first noting that 


f'V’L f/pripH t/erL 


wr 


(27) 


Then, using equations (18) gives 


£ 


K’W-Kn -K’\ 

\KT K'l Kl 


(28) 


and using a! 1 = A y 0) = ’ = A y 1 = 0 gives the desired result 


£ 


/p‘ M \ dz _ /o\ 

K7 W 


(29) 


Therefore, the constants G x and G are found by integrating equations (25) to get 



(30a) 
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With G x and G y known, equations (20) are written as 


F x '(z)=H x '(z) + ^-l 

C 55 

F/(z)=H/(z)+^-l 

r 

v — 44 


and constitutive equations (22) become 


1° y Z ( x ’ y> z ) \ 

\ O xz (x, y, z) f 


v^44 ^45 

pW pW 


H/(z) 

G 


Yy Z ( X ’ y) | 

+ — 

— (k) 

- 1 


C 44 


l 

H/(x) 

G x 

H — 

C (k) 

^55 

- 1 

Yl( x ’y) j 


(30b) 


(31a) 

(31b) 


(32) 


(It) 

The next step in the analysis is to obtain expressions for the zig-zag functions <F X (z) and 


4>y (z) in terms of the constants G x and G . Toward this goal, equations (18) are re-written as 


Al k ’ 

= A x kl) 


3'x k) 

A ^ 

A ( k - ! ) 

1 ( k )r 

,( k ) 

Ay 

= A y 

+ h (: 

E 


(33a) 

(33b) 


Examining these equations for successive values of the index k, starting with k = 1 and using 
a! 1 = A, 1 = 0 , it is seen by mathematical induction that 


. W V o (p) u 
A x =2j IE h 


(p) 


A (k) V Q (P) T, 

Ay — 2-k Py h (p) 

p= 1 


Substituting these expressions into equations (17) gives 


< hx k ’( z ) = 2 P ( x Pl h lp) + p'x k) (z - z^p) 

P = 1 


(33c) 

(33d) 


(34a) 
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(34b) 


= 2 Py P)h (P) + P'y M ( Z - Z (lc-.)) 

p = 1 

Now, consider h (k) = z (k) - z (k l) for successive values of the index k, starting with k = 1. Noting 
that z (<>) = - j , mathematical induction yields 


z ( k) - 9 + 


K 

2h 


(p) 


Thus, using equations (25) and (35), equations (34) are expressed as 


d>! k, (z) 


<’( z ) 


)\^-l 

k- 1 

+ g x 2 

u u 

^(p| — (k) 

|_ C55 

p= 1 

Lc 55 

C 55 J 

ll-^-il 

1 ^( k l 

k- 1 

+ G y 2 

U 3_1 

-^(P) — (k) 

/~1 

_ V-44 

P = 1 

r 

^44 

C 44 J 


* (p> 


* (p) 


(35) 


(36a) 

(36b) 


Let M denote the ply that corresponds to the plate midplane, z = 0. At z = 0, equations (36) 
reduce to 


*:>)=4 


G v 


C 


(M) 


- 1 


+ G X 2 


p = 1 


C 


(p) 


-|M| 


‘(P) 




+ G y 2 


1 




h 


(P) 


(37a) 

(37b) 


Examination of equations (14) indicates that the constraints F x (0) = F y (0) = 0 require 

H x (0) + or’(O) = 0 (38a) 

H v (o) + <>(o) = 0 (38b) 

Similarly, the constraints given by equations (12), which enforce the traction boundary conditions 
on the top and bottom plate surface, yield 


h.’H) 

h.'H) 



1 =0 
1 =0 


(39a) 

(39b) 
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(39c) 



Solving these equations yields 



Substituting these results into equations (40) gives 


(42a) 

(42b) 

(42c) 

(42d) 
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H,w = f,M - 4>I“’(0) ^ - f f Stt + f.'(- l)l( 1 - § ) - f f 3k * f.'(+ l)l( 1 + g ) (43a) 

C-5S C 55 


H,(z) = f,(z) - *:>) + z - § [5y + f/(- i)l( 1 - f ) - 1 f St + f/(+ l)l( 1 + f ) ( 43b > 

C44 C44 


In addition, 


H.'(z) = f,'w + 1 - ife + f.'(- |)l(i - f ) - !& + (.'(+ i)! 1 + t) <44a) 

^ 55 ^55 

H/(z) = f/(z) + 1 - I ^ + f/(- 1 ) (i - X ) " J Si + f y'( + l) f 1 + if ) (44b) 


In the present study, the functional form of f x (z) and f y (z) are obtained by requiring H x (z) and 
H y (z) to reduce to the distributions expected for a homogenous, orthotropic plate. For this very 

_(k) _(k) 

special class of plates, C 44 = G 23 and C 55 =G 13 , where G 13 and G 23 are the shear moduli. Thus, 
equations (30) give G x = G 13 and G y = G 23 . In addition, 4> X M) (0) = O y M| (o) = 0. Equations (43) 


H x (z)=f x (z) 

- ' 
2 x 

H) 

c-t: 

| + f x'(+|) 


H y (z) = f y (z) 

2 y 1 

(-1) 

(>-f) 

+ f/(+l) 

(■+§( 


and equations (44) reduce to 

H.'(z) = f.'(z) - |)(l - (f ) + f.'(+ l)( 1 + t) 


H,'(z) = f,'(z) - |)( 1 - f ) + f/( + |)(l + f ) 

Moreover, the transverse sheaing stresses given by equation (32) become 


(45a) 

(45b) 


(46a) 

(46b) 


/ cr yz (x,y,z) 
\ cr xz (x,y, z) 


G 23 0 j H/(z)y;(x,y) 
0 °13 | H/(z) y1(x, y) 


For a homogenous plate with traction-free shear-stress boundary conditions on the top and bottom 
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surfaces, the distributions of transverse-shear stresses are symmetric about the plate midplane. 
Thus, H/(- z) = H/(z) and H/(- z) = H/(z). Enforcing these symmetry requirements yields 


f» - f.'(- z) - f [(,'(+ |)-f.'(-|)]=0 (48a) 

f/M -f/t-z) - x[ f /(+ i) - f/(- 1)] = 0 (48b) 

for all values of -h/2 < z < h/2. Therefore, the functions f x (z) and f y (z) are required to satisfy: 

f x (0) = f y (0) = 0 (49a) 

fx'(-z) = f/(z) (49b) 

f y'(- z ) = f /( z ) (49c) 

for all values of -h/2 < z < h/2. Additional simplifications are obtained by specifying 

f.'RH'HH (so) 

Enforcing these conditions in equations (43) and (44) yields, without a loss in generality, 


and 



Using these expressions with equations (36), equations (14) are expressed as 


F x (z) = f x (z) + <'(z, 
F y (z)=f y (z) + <’(z) 


(51a) 

(51b) 


(52a) 

(52b) 
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(53a) 

(53b) 



where f x (z) and f (z) are continuous functions with continuous derivatives that satisfy equations 
(47) and (48). The functions ’J'Tfz) and 'I 1 ' (z) are the zig-zag enrichment functions given by 


^ k, (z) = <&‘ k, (z) - 4>r(0) + 

<(*) - <V) - O) + 




(54a) 

(54b) 


where < h! M (z) and ‘F '(z) are given by equations (36), ^ M, (o) and ^y (°) are given by 
equations (37), and G x and G y are given by equations (30). The derivatives of F x (z) and F y (z) 
are obtained by substituing equations (52) into equations (31); which yields, 


F/(z) = f/(z) + W, '(z 


F/(z) = f/(z) + ^’'(z 


where 



(55a) 

(55b) 


(56a) 

(56b) 


For the parabolic distribution of transverse shearing stresses commonly found in the technical 
literature for homogeneous plates, 


f *(z) = f y( Z ) = z [ 

f/(z)=f/(z) = 



These functions are illustrated in figure 9. 


(57a) 

(57b) 
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Examples 


Results obtained from equations (53) and (55), based on equations (57), are presented for two 
examples in this section. The first example is a [- 1 5/30/0/90] T laminated composite plate that 
exhibits the full extent of plate anisotropy. The second example is a sandwich plate with a 0.42- 
inch-thick isotropic core and identical [±45J laminated composite face plates, denoted by [±45 4 / 
core] s . For both examples, the principal transverse- shear moduli of a ply are given by G 13 = 0.864 
x 10 6 psi and G 23 = 0.368 x 10 6 psi, and the thickness of each ply is 0.005 in. For the core of the 
sandwich plate, G 13 = G 23 = G core = 0.100 x 10 6 psi. 

For the [- 1 5/30/0/90] T laminated composite plate, equations (30) yield G x /G 13 = 0.72 and 

G y /G 93 = 1.29. In addition, O '“'(()) and O,: M '(0) defined by equations (37) are given by -0.0021 
and 0.0008 in., respectively, where M = 2 and 3 contain the plate midplane and give identical 
results. Graphs of F x (z) and F y (z) given by equation (53a) and (53b) are shown in figures 10 and 
11, respectively. These graphs clearly show the effects of laminate asymmetry and inplane- 
coordinate direction on the inplane-displacement distributions. Moreover, both graphs depict 
piecewise-cubic curves that are continuous at the ply interfaces. Similarly, graphs of the 
derivatives of F x (z) and F y (z) given by equation (55a) and (55b) are also shown in figures 10 and 
11, respectively. These two derivative graphs depict piecewise-quadratic curves that are 
discontinuous at the ply interfaces, consistent with the refined zig-zag theory of Tessler et. al. 43 

For the seventeen-ply [±45 4 /core] s sandwich plate, the core comprises 84% of the plate 

thickness. Equations (30) yield G x /G core = G y /G core = 1.16. In addition, %"(()) and 4>' (o) 

defined by equations (37) are equal to zero, where M = 9 corresponds to the plate midplane. For 
this plate, F x (z) and F y (z) given by equation (53a) and (53b) are identical, and a graph of F x (z) is 
shown in figure 12. This graph which clearly shows the presence of laminate symmetry is 
composed piecewise-cubic curves that are also continuous at the ply interfaces, with the most 
variation occurring in the more flexible core, as expected. Note that no variation in F x (z) is shown 

_(k) _(k) — (k) 

across the faceplates, consistent with the fact that the , C 45 , and C 55 are identical for two 
adjacent +45 and - 45 degree plies. The corresponding graph of the derivatives of F x (z) is also 
shown in figure 12. This derivative graph also depicts piecewise-quadratic curves that are 
discontinuous at the ply interfaces. Moreover, the distribution of the curves is symmetric across 
the plate thickness. 


CONCLUDING REMARKS 

A detailed anaylsis has been presented that shows how to enrich the kinematics of classical 
Kirchhoff plate theory by using a set of continuous piecewise-cubic functions referred to in the 
technical literature as zig-zag functions. This work was motivated by the desire to obtain realistic 
estimates of the effects of laminate heterogeneity and asymmetry on the variations of the inplane 
displacements and transverse shearing stresses, for use with a {3, 0} plate theory in which these 
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distributions are specified apriori. The enrichment is based on the improved zig-zag functions, 
and the corresponding analysis methodology, presented recently by Tessler, Di Scuva, and 
Gherlone. With the approach that has been presented herein, the inplane displacements are 
represented by a set of continuous piecewise-cubic functions. In contrast, the transverse shearing 
stresses and strains are represented by a set of piecewise-quadratic curves that are discontinuous 
at the ply interfaces. Corresponding results have also been presented for a general asymmetric 
laminate and a sandwich plate with identical laminated composite face plates and an isotropic 
core. These results demonstrate clearly that the zig-zag enrichment captures the expected effects 
of laminate heterogeneity and asymmetry. 
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Fig. 2 Principal material coordinate system at z = constant. 
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Fig. 3 Kinematics of first-order shear-deformation plate theory in the x-z plane. 



Fig. 4 Deformation of a differential element in the x-z plane, based 
on first-order shear-deformation plate theory. 


22 


z 

n 


Undeformed 

configuration 




► x 


Plane y = constant 


Fig. 5 Deformation of plate cross-section, based on the kinematics used in the 
present study. 
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(a) Differential element at z = 0 (b) Differential element at z = constant 

Fig. 7 Deformation of a differential element in the x - z plane, based 
on the kinematics of the present study. 
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Fig. 9 Graphs of the function f x (z) and its derivative. 
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Fig. 12 Graphs of the function F x (z) and its derivative for [±45 4 /core] s sandwich plate. 
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